Higher-order multipole amplitude measurement in — )- 7Xc2 
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Abstract 

Using 106 X 10*^ V' events collected with the BESIII detector at the BEPCII storage ring, the 
higher-order multipole amplitudes in the radiative transition ■0' — > 7Xc2 'yTr^Tr~ /'yK^K~ are 
measured. A fit to the Xc2 production and decay angular distributions yields M2 = 0.046 ± 0.010 ± 
0.013 and £^3 = 0.015ib0.008ib0.018, where the first errors are statistical and the second systematic. 
Here M2 denotes the normalized magnetic quadrupole amplitude and E3 the normalized electric 
octupole amplitude. This measurement shows evidence for the existence of the M2 signal with 
4.4(7 statistical significance and is consistent with the charm quark having no anomalous magnetic 
moment. 

PACS numbers: 13.20.Gd, 13.25.Gv, 13.40.Hq 
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I. INTRODUCTION 



The radiative transitions ip' — )■ 7Xcj {J = 1? 2) provide information about the electromag- 
netic interaction between charm and anti-charm quarks in charmonia and allow investigation 
of many interesting topics, including whether the charm quark has an anomalous magnetic 
moment [l|, ^ or if there is S-wave and D-wave state mixing jsf. In general, the transi- 
tion amplitude of radiative decays of charmonium states is dominated by the electric dipole 
(El) contribution, with higher multipoles suppressed by powers of photon energy divided by 
charm quark mass [J] . The search for contributions of higher-order multipole amplitude is 
of interest as a source of information on the charm quark's magnetic moment; the possibility 
of anomalous magnetic moments of heavy quarks being larger than those of lighter ones was 
raised in Ref. j5|. In — )■ 7Xc2; taking the charm quark to have a mass = 1.5 GeV/c^ 
and an anomalous magnetic moment k, M2 = 0.029(1 + k) is predicted 0]. 

Disagreement between a pure El calculation and experimental measurements |6| hint 
that higher order multipole amplitudes may exist. These would be reflected in the angular 
distributions of both the radiative photon and the final state particles 0, Isf. Thus, careful 
investigation of the angular distributions is important. 

Several experiments, including the Crystal Ball experiment in if)' — )■ 7Xci,c2 ~^ 11 J ~^ 
77£"'"£~(£ = e or n) jof, the E-835 experiment in pp — >■ Xci,c2 — >■ iJ/i' — ^ 7e+e~ (lol |. 



the E-760 experiment in pp — >_yc2 — ^ iJ/i^ 7e+e~ [ll|, and the BESII experiment in 
ip' — )• 7Xc2 jTi^Ti' /'jK^ K~ [l2| have searched for higher-order multipole amplitudes. Due 
to their limited statistics, they were unable to provide evidence for the existence of higher- 
order multipoles. More recently, the CLEO experiment reported measurements of higher- 



order multipole amplitudes in ip' — )■ iXci,c2 — ^ llJ/i^ 77£^£~ [13|, where significant 
M2 contributions were found in the ip' — )■ iXci and Xci,c2 — ^ iJ/i^ transitions. Tables U 
and HI] summarize the experimental measurements on searches for higher-order multipole 
amplitudes. 

TABLE I: Current experimental measurements of the normalized M2 contributions in the decays 
Xci -^iJl^ and V' lXc\- 



Experiment 


Xci iJ/i^ 




IXci Sij 


^nal Events 


Crystal Ball [9] 


U.UUZ_Q g2Q 


0.077 


+0.050 
-0.045 


921 


E-835 [10] 


0.002 ±0.032 ±0.004 






2090 


CLEO-c ^ 


-0.0626 ± 0.0063 ± 0.0024 0. 


0276 ± 0.0073 ± 0.0023 


39363 



In this article, (1.06 ± 0.04) x 10^ ip' events IJ] accumulated in the BESIII experiment 
are used in the selection of ip' — 7Xc2,Xc2 — ^ 71^11' /K^K~ events, which allow the deter- 
mination of the higher-order multipole amplitudes in the ip' — )■ 7Xc2 transition. 



II. THE BESIII EXPERIMENT AND DATA SET 

This analysis is based on a 156.4 pb^^ of ip' data corresponding to (1.06 ± 0.04) x 10^ ip' 



events [IJ] collected with the BESIII detector [15| operating at the BEPCII Collider [16 
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TABLE II: Current experimental measurements of the normalized M2 contributions in the decays 
Xc2 -^^J/Tp and V' 7Xc2- 



Experiment 


Xc2 ^iJ/i} 


V'' ^ 7Xc2 Sij 


^nal Events 


Crystal Ball [9] 


n oqq+U.116 
U.OOO_Q 


n 1 q9+U.U98 

U.iOZ_Q 


441 


E-760 [11] 


-0.14 ±0.06 




1904 


E-835 [10] 
BESII [12] 


-0.093l[]:[}^? ± 0.006 


U.U01_o 036 


5908 
731 


CLEO-c [13] 


-0.079 ± 0.019 ± 0.003 0.002 ± 0.014 ± 0.004 


19755 



In addition, an off-resonance sample of 42.6 pb ^ taken at ^/s = 3.65 GeV is used for the 
study of continuum backgrounds. 

BESIII/BEPCII [lH is a major upgrade of the BESII experiment at the BEPC accel- 
erator [l6[ for studies of hadron spectroscopy and r-charm physics [17]. The design peak 
luminosity of the double-ring e"'"e~ collider, BEPCII, is 10^^ cm~^s~^ at beam currents of 
0.93 A. The BESIII detector with a geometrical acceptance of 93% of An, consists of the 
following main components: 1) a small-celled, helium-based main draft chamber (MDC) 
with 43 layers. The average single wire resolution is 135 /^m, and the momentum resolu- 
tion for 1 GeV/c charged particles in a 1 T magnetic field is 0.5%; 2) an electromagnetic 
calorimeter (EMC) made of 6240 Csl (Tl) crystals arranged in a cylindrical shape (barrel) 
plus two end-caps. For 1.0 GeV photons, the energy resolution is 2.5% in the barrel and 5% 
in the endcaps, and the position resolution is 6 mm in the barrel and 9 mm in the end-caps; 
3) a Time-Of-Flight system (TOP) for particle identification (PID) composed of a barrel 
part made of two layers with 88 pieces of 5 cm thick, 2.4 m long plastic scintillators in each 
layer, and two end-caps with 48 fan-shaped, 5 cm thick, plastic scintillators in each end-cap. 
The time resolution is 80 ps in the barrel, and 110 ps in the end-caps, corresponding to a 
K/tt separation by more than 2a for momenta below about 1 GeV/c; 4) a muon chamber 
system (MUG) made of 1000 m^ of Resistive Plate Chambers (RPG) arranged in 9 layers 
in the barrel and 8 layers in the end-caps and incorporated in the return iron yoke of the 
superconducting magnet. The position resolution is about 2 cm. 

The optimization of the event selection and the estimation of physics backgrounds are 
performed through Monte Garlo (MG) simulations. The GEANT4-based simulation soft- 



ware BOOST [18| includes the geometric and material description of the BESIII detectors 
and the detector response and digitization models, as well as the tracking of the detector 
running conditions and performance. The production of the ip' resonance is simulated by 



the MG event generator KKMG 19(], while the decays are generated by EVTGEN |2(| for 
known decay modes with branching ratios being set to PDG [2l| world average values, and 
by LUNDGHARM [22| for the remaining unknown decays. The analysis is performed in 
the framework of the BESIII offline software system |22] which takes care of the detector 
calibration, event reconstruction and data storage. 

MG samples of ip' — ^ 7Xco,c2 — ^ ■yn~^7i~ /'jK~^K~ are generated according to phase space to 
determine the normalization factors in the partial wave analysis [l2[|23l, and MG samples of 
ip' — )■ (7)e+e~/(7)/i+/i~ and ip' ~^ XJ/ip (X = vr^vr^,?]) with J/ip — )■ {j) fj,~ are generated 
for background studies. 
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III. DATA ANALYSIS 



Charged tracks are reconstructed in the MDC, and the number of charged tracks is 
required to be two and have no net charge. For each track, the polar angle must satisfy 
I cos6'| < 0.93, and the point of closest approach must be within ±10 cm of the interaction 
point in the beam direction and within ±1 cm of the beam line in the plane perpendicular 
to the beam. The TOF (both end-cap and Barrel) and dE/dx measurements for each 
charged track are used to calculate Xp/d(0 values and the corresponding confidence levels 
Probpi£,{i) for the hypotheses that a track is a pion, kaon, or proton, where i {i = n/K/p) 
is the particle type. For pion candidates, ProbpjD^n) > 0.001 is required, while for kaon 
candidates, ProbpjD^K) > ProhpioiT^) and Probpj£){K) > 0.001 are required. 

Electromagnetic showers are reconstructed by clustering EMC crystal energies. The 
energy deposited in nearby TOF counters is included to improve the reconstruction effi- 
ciency and energy resolution. Showers identified as photon candidates must satisfy fidu- 
cial and shower-quality requirements. The minimum energy is 25 MeV for barrel showers 
(I cos^l < 0.8) and 50 MeV for end-cap showers (0.86 < | cos^^| < 0.92). The showers in the 
angular range between the barrel and end-cap are poorly reconstructed and excluded from 
the analysis. To eliminate showers from charged particles, a photon must be separated by 
at least 20 degrees from any charged track. EMC cluster timing requirements are used to 
suppress electronic noise and energy deposits unrelated to the event. The number of good 
photon candidates is required to be larger than or equal to one in each event, and the photon 
with the highest energy is regarded as the radiative photon from ip' — )• 7Xc2- 

In order to separate pions and kaons more effectively and to distinguish Xco and Xc2 more 
clearly, a four-constraint kinematic fit (4-C fit) is performed with the two charged tracks and 
the radiative photon candidate under the hypotheses that the two tracks are either tt+tt" 
or K^K~, and the kinematic chi-squares, xl and x'k^ are determined. If xl < Xk 
xl < 60, the event is categorized as 77r"''7r~; otherwise, if Xk < xi and Xk < 60, the event 
is categorized as 'yK^K~. For the selected 77? "•"vr" and 'yK^K'' candidate events, at least 
one of the charged tracks is required to be identified as a vr for 77r+7r~ or a K for '~fK^K~. 

To remove e+e^ {■y)e^e~ and ip' {■y)e~^e~' backgrounds, the deposited energy of 
each track in the EMC is required to be less than 1.4 GeV, and the observed ionization 
is also required to be within 3a of the expected value for each track. Furthermore, in 
77r"'"7r~, observed ionization of charged tracks is required to be within 2a of the expected 
value when the polar angle of the charged track is within the EMC insensitive region (i.e. 
0.81 < I cos^l < 0.86). These requirements remove almost all events with two electron tracks 
but still keep the efficiencies for the signal channels very high: 96% for 77r+7r^ and 97% for 

In 77r+7r~, there are e~^e~ — )■ {j) fj,'^ fj,~ and ip' ~^ {1)1^^ backgrounds due to n/ix 
misidentification. In order to remove the backgrounds with the deposited energy in 

the EMC of at least one of the charged tracks is required to be larger than 0.34 GeV. Over 
99% of the e+e^ — )■ (7)/!+/!^ and ip' — )■ (7)/i+yU~ backgrounds are removed after applying 
this requirement. Since ^/K misidentification is quite small, it is not necessary to apply 
this requirement in the jK^K' decay. 

After performing all the above selection criteria, clean ip' — )■ 7Xco,c2 — ^ jTr'^ir' /'yK~^K~ 
data samples are obtained. The tt+tt" and K^K~ invariant mass distributions are shown 
in Fig. [TJ Clear Xco and Xc2 signals are observed, and the background levels within the Xc2 
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FIG. 1: The invariant mass distributions of vr^vr" (left) and K^K^ (right) for the selected "fTi~^iT~ 
and 'yK~^K~ events from ip' data. Dots with error bars are data while blank histograms are the 
sum of MC simulated backgrounds and normalized continuum background (estimated from the 
data sample taken at ^/s = 3.65 GeV). 



signal region between 3.53 and 3.59 GeV/c^ are 2.6% [0.7% ip' (7)/i+/i", 0.8% normalized 
continuum, 0.3% cross contamination from Xc2 K^K~, 0.7% Xco tciil, and 0.1% ip' — )■ 
7r+7r^/7r+7r^7r° events] for 77r+7r^ and 2.1% [0.7% cross contamination from Xc2 — >■ tt+tt^, 
1.1% Xco tail, and 0.3% ip' — )■ K^K~ events] for 'yK^K^ . The highest mass peak corresponds 
to ip' decays to two charged tracks that are kinematically fitted with an unassociated low 
energy photon. Requiring the invariant mass of the two charged tracks to be between 3.53 
and 3.59 GeV/c^ to select Xc2, 7154 ■j7i~^7i~ events and 6657 jK'^K' events are obtained. 



IV. FIT TO THE ANGULAR DISTRIBUTIONS 



The formulae for the helicity amplitudes in ip' — )■ 7Xc2 — ^ 7PP (P =_[[/ K) , which include 
higher-order multipole amphtudes, are widely discussed in Refs. 

3 

+ -(1 + cos2^^)sin^^M2/^ 



V2 



sin(26'^) sin(26'M) sin^ 9m cos (pMxy 



+ a/3 sin(26'-y) sin(26'M)(3 cos^ 6m — 1) cos 0^3^ 
+ sin^ sin^ 6'm(3 cos^ Qm — 1) cos 20^/1/ 
+ (1 + cos2^T,)(3cos2eAf - 1)^ 



(1) 



where x = Ai/Aq, y = A2/A0, and v4o,i,2 are the Xc2 helicity 0, 1, 2 amplitudes, respectively. 
6^ is the polar angle of the radiative photon where the electron beam is defined as the z axis 
in the e^e~ center-of-mass frame, and 6m and (pM are the polar and azimuthal angles of 
the n/K in the Xc2 rest frame, where the polar axis is defined with respect to the radiative 
photon direction and 0m = is defined by the electron beam direction. 
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An unbinned maximum likelihood fit to the joint production and decay angular distribu- 
tion is performed to determine x and y values. We define six factors [24ii|: 



ai = 3sm'^9^sm\29M), (2) 
3 

02 = - (1 + cos^^^) sin^ 6'm, (3) 
3 

as = ^sin(2^^)sin(26'M)sin^6'MCOs0Af, (4) 

a4 = \/3 sin(26'-y) sin(2^^M)(3 cos^ — 1) cos 0M, (5) 

as = a/6 sin^ 6^^ sin^ 6'm(3 cos^ — 1) cos20jv/5 (6) 

as = (1 + cos2^^)(3cos2^M - 1)^. (7) 



The mean values of ai, a2, 03, 04, as, ag are determined with ip' — )■ 7Xc2,Xc2 — ^ PP MC 
events, where phase space is used for the simulation of all the angular distributions: 

an = — , n = l,...,6 (8) 

Here N is the number of events after all selections from phase space MC samples. We 
integrate firstly the parts independent of the parameters in the angular distribution to make 
the fit faster. Since an is calculated with phase space MC events, thus it accounts for the 
detector acceptance effects naturally. 

Then, the constructed probability-density function (pdf) is written as: 

aix^ + a2|/^ -l- a^xy + a^x + a^y + 

In practice, we use the log-likelihood function, which is given by ln£ = S^^^ In for 
convenience, where the sum is over the events in the signal region. The dominant background 
events are simulated by MC events and their normalized contributions are subtracted in In C 
value, i.e. ln£s = ln£ — ln£{,, where \nCh is the normalized sum of logarithmic likelihood 
values from background events. 

Before fitting to the data, input and output checks have been done using MC samples, 
and the checked results verify the validity of the fitting procedure. An unbinned maximum 
likelihood fit to the 77r"^7r~ and -^K^ production and decay angular distributions yield 

x^ = 1.55^°:°?, y^ = 2Mt^^Z Pn = 0.890. (10) 
xk = 1.55 ± 0.08, yK = ^.U^oH, Pk = 0.902. (11) 

where the errors are statistical, and p^, px are the correlation coefficients between x and y 
for 771"*" 7r~ and 'yK~^K~ , respectively. A simultaneous fit to 771 "•"vr" and '-)K^K~ gives 

X = 1.55 ± 0.05, y = 2.10 ± 0.07, p = 0.896 (12) 

where errors are statistical, and p is the correlation coefficient between x and y. The nor- 
malized M2 and E3 amplitudes are calculated to be [l|: 

M2 = 0.046 ± 0.010, E3 = 0.015 ± 0.008 (13) 
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based on the simultaneous fit result, where errors are statistical only. Figure |2] shows the 
angular distributions of data and the fitted results both for 77?"*" vr" and 'yK'^K'^ events. 
Good agreement is observed for all angular distributions. 

The goodness of the fit is estimated using Pearson's t^st 25|. The data sample is 
divided into 8x8x8 = 512 bins in cos6'-y, cos^m, and (f)M, and the value is calculated 
as: 



,2 



pT_„MC; 



E '"' :Ot' ' (14) 



where nf^ and nf^*^ are the number of the observed events in i-th bin from data and the 
corresponding number of normalized events from MC using x and y fixed to the values 
determined in the analysis. Here MC events are 20 times more than data events. For 
bins with less than 7 events, we add the events into the adjacent bin. The result yields 
Xl/n.d.f = 377.5/368 = 1.03 for -fn+n- and x\/n.d.f = 348.1/354 = 0.98 for -fK+R-, 
where n.d.f is the number of the degrees of freedom. These results show that the fits are 
good. 



V. SYSTEMATIC ERRORS 



A. MC simulation of detector response 



The consistency between data and MC simulation for Xc2 events can be tested using Xco 
events. The angular distribution of Xco is unambiguous, i.e., Wq = 1 + cos^ 9^. If we replace 
the (3 cos^ — 1)^ term in Eq. 1 by 1, then Eq. 1 becomes equal to Wq when both x and y 
equal zero. Therefore, if we fit the angular distribution of Xco events with a modified Eq. 1 
using the same method as in Xc2 decays, x = and y = are expected. Non-zero x and y 
values from the fit refiect the difference between data and MC and give a measure of the 
systematic error due to the MC simulation of the detector response. The fitted results are 
x^ = 0.049l|^-[!}?, y^ = -(0.024 ± 0.011), = 0.037 for -fn+n-, xk = 0.073 ± 0.015, vk = 
0.004 ± 0.010, pk = 0.077 for jK+R-, and x = 0.062 ± 0.011, y = -(0.008 ± 0.007), p = 
0.058 for the simultaneous fit. The systematic error is taken as the shift from plus its error. 
Assuming the correlation factor is 1 between x and y for central value, 0.06, 0.04, and —0.63 
are obtained for Ax,r, ^yn and p^^"*, respectively; 0.08, 0.02, and 0.52 are obtained for Axk, 
Ayx and p^'^, respectively; and for the simultaneous fit, 0.07, 0.02, and —0.43 are obtained 
for Ax, Ay and p*^'^, respectively. Studies with the MC demonstrate that a systematic error 
in modeling the 9^, (j)M, or 9m efficiency produces a shift of x and y of approximately the 
same size in both Xco and Xc2 samples, where the latter sample is generated with our nominal 
results for x and y. Therefore, we assume the observed shift from x = and ?/ = for the 
true Xco data is an estimate of the systematic error on the measured values of x and y for 
radiative decays to Xc2- 

The systematic error of the detector response contains systematic errors associated with 
the simulation of charged track finding, photon detection efficiency, mass resolution of Xc2, 
kinematic fit, PID efficiency, trigger efficiency, etc. Comparisons between data and MC 
simulation for the angular distributions of Xco decay events are shown in Fig. [3l 
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FIG. 2: The angular distributions of cos9^, cosOm and 4>m for Xc2 — ^ 7r'''vr (left) and Xc2 
K^K~ (right), where the dots with error bars are data and the histograms are the fitted results. 



B. ^J' ^J' background in 777+77 



In 77r"^7r~, the dominant backgrounds are {'~f)fi^fi~ events and continuum events, which 
contribute at the 1.5% level. The ip' — {'~f)fi^fi" background events are estimated by 
MC simulation, while the continuum events are estimated using the data sample taken at 
^/s = 3.65 GeV. The /x/tt misidentification ratio has been checked using a control sample of 
ip' —7- 71^71^ J/'ip —7- TT+TT^/i+yU^ and its value is measured to be (1.7±2.6)%. Considering the 



uncertainty from tp' — fi~^fi branching raito (10.4%) 2l|], the {'j)fi'^fi background level is 



determined to be (1.5 ± 0.2)%. We vary the background level by la (from 1.5% to 1.7% or 
1.3%) in the fit and take the difference of the fitted x and y values as the systematic error 
due to the /i+/i~ background uncertainty. The differences are Ax = 0.02, Ay = 0.03. 
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FIG. 3: Comparisons of cos 9-y, cos 9m , and (pM angular distributions between data (dots with error 
bars) and MC simulation (histograms) for XcO — ^ tt+tt" (left) and XcO — ^ K^K~ (right). 



C. 777"*" TT and 'yK^ K cross contamination 



The systematic error arising from 77r+7r~ and jK^K' cross contamination is also deter- 
mined by MC simulation. The 'yK~^K~ background contamination in 771 "•"vr" is 0.3%, while 
the 77r+7r~ background contamination in 'yK~^K~ is 0.7%. Signal MC samples of 771 "•"tt" 
and '^K^K~ with x = a/3, y = a/6 are generated and mixed according to the estimated 
amount of cross contamination determined by MC simulation. The differences on the fitted 
X and y values are taken as the systematic errors due to the 77r+7r^ and '-yK^K~ cross con- 
tamination, which are Ax^ = 0.01, ^y-^ = 0.02 for 77r"'"7r~ and Axk = 0.01, Ayx = 0.02 
for -fK+K-. 
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TABLE III: Summary of the systematic errors and correlations. 



Source x.^ xk Vk Pk 

MC simulation 0.06 0.04 -0.63 0.08 0.02 0.52 

background 0.02 0.03 1 - - - 

tt/K cross contamination 0.01 0.02 1 0.01 0.02 1 
XcO contamination . _ _ Q.Ol 0.02 1 

Total 0.07 0.06 -0.17 0.08 0.04 0.38 



D- XcO contamination 



There are some Xco events in the Xc2 signal region due to the overlap of the tt^ti' / K^K~ 
invariant mass peaks. The contamination from Xco t^^t^~ background events is only about 
0.7% for Xc2 TT^TT~ , and the contamination from Xco K^K~ background events is 1.1% 
for Xc2 K^K^ according to MC simulations. Signal MC samples of Xco ~^ ti^ti^ / K^K~ 
are generated and mixed into Xc2 ti^ti^ / signal MC samples according to the 
estimated contamination ratio determined by MC simulation. The differences on the fitted 
X and y values from the input ones are taken as the systematic errors. The systematic 
uncertainties due to Xco contamination can be ignored for Xc2 ir^n', while they are 
Axk = 0.01, AyK = 0.02 for Xc2 ^ K+R-. 



E. Total systematic error 

The systematic error sources discussed above are summarized in Table IIIII Here the cor- 
relation coefficients (p^ and pk) from background uncertainties including fi^fi~ background, 
71 /K cross contamination background, and Xco background contamination are set to be 1, 
and the total correlation coefficient p is calculated as p = Ei^jej—mi where i runs over 

all systematic errors. The total systematic errors are Ax^ = 0.07, Ay^ = 0.06 in 77r"'"7r~, 
Axx = 0.08, AyK = 0.04 in -fK+K', and Ax = 0.07, Ay = 0.05 in the simultaneous fit. 



VI. CONCLUSION AND DISCUSSION 

The final helicity amplitude results are: 

x^ = 1.55l°:°? ± 0.07, y^ = 2.06^°;^° ± 0.06, pf"* = 0.890, p^J^' = -0.17 (15) 
for 77r+7r^, 

Xk = 1.55 ± 0.08 ± 0.08, yK = ^.U^oH ± 0.04, p^'^* = 0.902, p^" = 0.38 (16) 

for '-fK^K^ , and 

X = 1.55 ± 0.05 ± 0.07, y = 2.10 ± 0.07 ± 0.05, p"*"* = 0.896, p'^' = 0.26 (17) 
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FIG. 4: Experimental measurements for the normalized M2 amplitude together with theoretical 
prediction [l], [3] assuming the charm quark mass to be 1.5 GeV/c^ and no anomalous magnetic 
moment. CLEO's result is from free E3 amplitude fit. 



from the simultaneous fit, where the first errors are statistical and the second systematic. 
pstat psys ^YiQ correlation coefficients between x and y of the statistical and systematic 
errors. Then the normalized M2 and E3 amplitude are determined to be: 

M2 = 0.046 ± 0.010 ± 0.013, E3 = 0.015 ± 0.008 ± 0.018 (18) 

where the first errors are statistical and the second systematic. By investigating the differ- 
ence of the logarithmic likelihoods between a pure El transition and the best nominal fit, 
the statistical significance of the M2 amplitude contribution is estimated to be 4.4(T, which 
means evidence of the M2 contribution has been observed. As for the E3 signal, the current 
measurement is consistent with zero. The M2 experimental results from different measure- 
ments are shown in Fig. HI Our measurement agrees with prediction and is consistent with 
CLEO's result within 2(t when E3 is free [3 and BESII's result |l2| within 1.7(t. 

In summary, the higher order multipole amplitudes in ip' — )■ 7Xc2 ^iti^ti^ /^K^K~ are 
studied with the BESllI experiment based on (1.06 ± 0.04) x 10^ ip' events. Evidence of an 
M2 amplitude is observed. This measurement agrees with prediction and is consistent with 
the charm quark having no anomalous magnetic moment [J, 0, 0] • 



VII. ACKNOWLEDGEMENT 



The BESIII collaboration thanks the staff of BEPCII and the computing center for their 
hard efforts. This work is supported in part by the Ministry of Science and Technology of 
China under Contract No. 2009CB825200; National Natural Science Foundation of China 
(NSFC) under Contracts Nos. 10625524, 10821063, 10825524, 10835001, 10935007; the 
Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; CAS under 
Contracts Nos. KJCX2-YW-N29, KJCX2-YW-N45; 100 Talents Program of CAS; Istituto 
Nazionale di Fisica Nucleare, Italy; Siberian Branch of Russian Academy of Science, joint 



13 



project No 32 with CAS; U. S. Department of Energy under Contracts Nos. DE-FG02- 
04ER41291, DE-FG02-91ER40682, DE-FG02-94ER40823; University of Groningen (RuG) 
and the Helmholtzzentrum fuer Schwerionenforschung GmbH (GSI), Darmstadt; WCU Pro- 
gram of National Research Foundation of Korea under Contract No. R32-2008-000-10155-0 



[1] J. L. Rosner, Phys. Rev. D 78, 114011 (2008). 

[2] K. J. Sebastian, H. Groteh and F. L. Ridener, Phys. Rev. D 45, 3163 (1992). 
[3] P. Moxhay and J. L. Rosner, Phys. Rev. D 28, 1132 (1983). 

[4] G. Karl, S. Meshkov, and J. L. Rosner, Phys. Rev. Lett. 45, 215 (1980); V. A. Novikov et al, 

Phys. Rept. 41, 1 (1978). 
[5] D. A. Geffen and W. Wilson, Phys. Rev. Lett. 44, 370 (1980). 
[6] E. Eichten et al, Phys. Rev. D 21, 203 (1980). 
[7] P. K. Kabir and A. J. G. Hey, Phys. Rev. D 13, 3161 (1976). 
[8] C. Edwards et al, Phys. Rev. D 25, 3065 (1982). 

[9] M. Oreglia et al (Crystal Bah Collaboration), Phys. Rev. D 25, 2259 (1982). 
[10] M. Ambrogiani et al (E835 Collaboration), Phys. Rev. D 65, 052002 (2002). 
[11] T. Armstrong et al (E760 Collaboration), Phys. Rev. D 48, 3037 (1993). 
[12] M. Abhkim et al (EES Cohaboration), Phys. Rev. D 70, 092004 (2004). 
[13] M. Artuso et al (CLEO Collaboration), Phys. Rev. D 80, 112003 (2009). 
[14] M. Ablikim et al (BESIII Collaboration), Phys. Rev. D 81, 052005 (2010). 
[15] M. Ablikim et al (BESIII Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 614, 
345 (2010). 

[16] J. Z. Bai et al (BES Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 344, 319 

(1994); 458, 627(2001). 
[17] D. M. Asner et al, Int. J. Mod. Phys. A 24, No. 1, 499 (2009). 

[18] Z. Y. Deng et al, HEP & NP 30, 371 (2006). 

[19] S. Jadach, B. F. L. Ward and Z. Was, Comp. Phys. Commu. 130, 260 (2000); S. Jadach, B. 

F. L. Ward and Z. Was, Phys. Rev. D 63, 113009 (2001). 
[20] R. G. Ping et al, HEP & NP 32, 599 (2008). 

[21] K. Nakamura et al (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[22] J. C. Chen, G. S. Huang, X. R. Qi, D. H. Zhang, and Y. S. Zhu, Phys. Rev. D 62, 034003 

(2000). 

[23] W. D. Li, H. M. Liu et al, in proceeding of CHEP06, Mumbai, India, 2006 edited by Sunanda 

Banerjee (Tata Institute of Fundamental Research , Mumbai, 2006). 
[24] Shi Xin, Master thesis, Peking University (2004). 

[25] W. Eadie et al, Statistical methods in experimental physics, Amsterdam-London: North- 
Holland Publishing Company, 1971. 



14 



